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When time-reversal symmetry is broken, quantum coherent systems with and without spin rota- 
tional symmetry exhibit the same universal behavior in their electric transport properties. We show 
that spin transport discriminates between these two cases. In systems with large charge conduc- 
tance, spin transport is essentially insensitive to the breaking of time-reversal symmetry, while in 
the opposite limit of a single exit transport channel, spin currents vanish identically in the presence 
of time-reversal symmetry but can be turned on by breaking it with an orbital magnetic field. 
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Introduction. Fifty years ago, Dyson showed that 
ensembles of unitary matrices that are invariant under 
general symmetry groups reduce to the direct product 
of three irreducible ensembles These three circular 
ensembles are labelled by an index (3 = 1,2,4 and are 
respectively invariant under the transformations 



S — > U T SU , orthogonal ensemble, (3 
S —> USV , unitary ensemble, (3 = 2, 



1. 



(la) 
(lb) 



S -> W R SW , symplectic ensemble, (3 = 4, (lc) 

where S is an element of the ensemble, U and V are 
arbitrary unitary matrices, W is a quaternion unitary 
matrix, U T is the transpose of U and W R = cr^'Wcr^' 
is the dual of W @- Here and below, cr^ , fi = x, y, z is a 
Pauli matrix. This classification carries over to electronic 
quantum transport, where the three classes are defined by 
the symmetry of the system Q . Systems without time- 
reversal symmetry (TRS) have a scattering matrix in the 
(3 — 2 ensemble, while systems with TRS arc differenti- 
ated by the presence ((3 = 1) or the absence ((3 = 4) of 
spin rotational symmetry (SRS). When TRS is broken, 
breaking SRS doubles the size of the scattering matrix, 
but does not generate a new ensemble. 

Quantum corrections to electric transport depend on 
the symmetry index f3, but are independent of the size 
N of the scattering matrix (giving the total number of 
transport channels from and to the scatterer) for large 
N 0. Recent investigations of spin transport showed 
that the spin conductance 



(2) 



constructed from the transmission block Sij of the scat- 
tering matrix connecting terminals i and j, also ex- 



hibit a character of universality 
Wi(Ni - l)Nj/N{2N - 1){2N 
Ni 

system and terminals i,j, and N 



ij 

Here, 

z j gives the number of transport channels between the 



in that varT, 



3) for /3 



Ei^i. The spin 



(a0\ 



= 



conductance fluctuates about zero average, (7 i 
and the resulting, typically nonzero spin current is gen- 
erated by the presence of the SRS breaking field. In the 



f3 = 4 ensemble one usually takes the latter field as spin- 
orbit interaction (SOI). In the absence of SOI, one has 
= 0. This is the case for (3 = 1 and, if Dyson's three- 
fold way applies to spin transport, for (3 = 2. In this 
manuscript we demonstrate that spin transport discrimi- 
nates between systems with and without SRS even when 
TRS is broken. Accordingly, a novel kind of universal- 
ity emerges in systems with broken SRS and TRS, with 
charge transport properties given by those of the (3 = 2 
ensemble, but with specific spin transport properties. 
The latter are similar to those of the (3 = 4 ensemble at 
large N, a finding already reported in Ref. Q for specific 
four-terminal setups, but deviate from it at small N. Our 
finding does not invalidate Dyson's classification- the lat- 
ter gives a complete classification of unitary scattering 
matrices and unless one introduces chiral or particle-hole 
symmetries [1, there is no new ensemble to be found. 
Instead our point is that spin-dependent obscrvablcs de- 
fine two sub-ensembles of the (3 = 2 ensemble, depending 
on whether they commute or not with the scattering ma- 
trix. In other words, we find that while universality in 
charge transport is affected only by the antiunitary sym- 
metries, universality in spin transport depends on both 
antiunitary and unitary symmetries. 

The model. We consider a mesoscopic conductor 
connected to any number of external electron reservoirs. 
There is no ferromagnetic exchange anywhere in the sys- 
tem, nor is there spin accumulation in the reservoirs and 
we neglect spin relaxation in the terminals. The mag- 
netoelectrically generated spin current due to the pres- 
ence of SOI inside the cavity is determined by the spin- 
dependent transmission coefficients of Eq. ([2]). For in- 
stance, in the simple case of a two-terminal setup, the 
generated spin current in the right lead along the polar- 
ization axis /j, = x, y, z is given by 



(3) 



with the voltage bias V applied across the sample. 

Semiclassical calculation. We first calculate the av- 
erage and mesoscopic fluctuations of the spin transmis- 
sion coefficients using the semiclassical theory of trans- 
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extended to take spin transport into ac- 
We write [U 



')Tr[U~<rMul] 



(4) 

The sums run over all trajectories starting at yo on a 
cross-section of the injection lead j and ending at y on 
the exit lead i. Trajectories have a stability given by Ay, 
which includes a prefactor (27rzTi) _1 / 2 as well as a Maslov 
index, and 5 7 gives the classical action accumulated on 
7, in units of h. SOI is incorporated in the matrices U 1 . 
The average spin conductance has been calculated semi- 
classically in Ref. 13( . In the absence of SOI, spins do not 



rotate, U 1 = <r<°> is the identity matrix, and one trivially 
obtains 7^ = 0. The leading-order approximation is to 
consider J7 7 G SU(2), where SOI rotate the spin of the 



electron along unperturbed classical trajectories [l2, 14 
In this manuscript, we will use this approximation be- 
cause, even thou gh i t neglects the geometric correlations 
reported in Ref. }13| . it is appropriate for our search of 
universality. At that level, the average spin conductance 
vanishes, (3y ) se micl = [3], which agrees with the ran- 
dom matrix theory (RMT) result of Ref. 

Having established that the average spin conductance 
vanishes regardless of the presence or absence of TRS and 
SRS, we next calculate spin conductance fluctuations. 
The leading-order diagrams contributing to var[T^Jscmici 
are shown in Fig. [TJ They are the same as those 
contributing to the (charge) transmission fluctuations 
[substituting <jM -> in Eq. ©]. In this case, 

Ref. [lH found that contributions c), d) and e) cancel 
out, furthermore, contribution b) vanishes upon break- 
ing of TRS. This can be achieved via a magnetic flux 
piercing the diagram's loop. From Fig. [TJ we see that 
contribution b) is the only one that is flux-sensitive, 
because the blue (dark) and the red (light) trajecto- 
ries accumulate the same flux-phase. From a semi- 
classical point of view, this is the origin of the halv- 
ing of the universal conductance fluctuations upon TRS 
breaking Extending this calculation to var[T^] somic i, 
we obtain that contributions a), b) and c) are multi- 
plied by a spin-dependent term TrlU^U^a^U^U^} x 
Tr[f/^ 3 J7^ 4 (T^^C/ 7l C/ 7l ], while contributions d) and e) are 
multiplied by |Tr[U'+(T^U 72 ] | 2 . All these terms vanish 
in the absence of SOI. In the presence of SOI, we evaluate 
them by averaging them over a uniform distribution of all 
CZy's over the SU(2) group, corresponding to totally bro- 
ken SRS. Following the standard procedure of performing 
orbital averages and spin averages separately, we obtain 
that, when SRS is totally broken, contributions a), 6) 
and c) acquire a prefactor ((•••)su(2) indicates an homo- 
geneous average over the SU(2) group) 

(Trl^t/t^^^^lTr^t/t^^^^Dsu^ =0, 

(5) 




Figure 1: (Color online) Semiclassical diagrams determining 
the conductance and spin conductance fluctuations to leading 
order in the number N ^> 1 of transport channels. Blue 
(dark) and red (light) trajectories travel in opposite direction 
in diagram b), which consequently vanishes in the presence of 
a large magnetic flux piercing the loop. All other diagrams 
are insensitive to the breaking of time-reversal symmetry. 

and thus vanish identically, while contributions d) and e) 
are multiplied by 

(|Tr[U-t i >)C/ T2 ]| 2 ) su(2) = l. (6) 

We conclude that the semiclassical contributions to the 
spin conductance fluctuations are those with a correlated 
encounter at the exit terminal, which in particular has 
the consequence that contribution b) vanishes. 

We obtain the variance of the spin conductance coeffi- 
cients as the sum of contributions d) and e), i.e. 



var^f 



scmicl 



NiNjN - TV, TV 2 /TV 3 



(7) 



The key point is that this result holds both in the ab- 
sence and in the presence of TRS, because both relevant 
contributions d) and e) arc sensitive neither to magnetic 
fluxes piercing their loops, nor to orbital magnetic field 
effects that do not alter the ergodicity of the classical 
trajectories. Thus, Eq. (|6j) gives the leading-order semi- 
classical expression for var G M for systems without SRS 
(with SOI) in both cases of conserved or broken TRS, 
as well as in the intermediate regime of partially broken 
TRS. Therefore, to leading order in the number TV 1 
of transport chanels, spin conductance fluctuations are 
insensitive to the breaking of TRS. In the next section, 
this result is confirmed using RMT. 

Random matrix theory calculation. We next use 
the method of Ref. [IH to calculate the RMT average 
and fluctuations of the spin conductance. We write [1] 
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otherwise , 



(8a) 
(8b) 

(8c) 
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where m and n are channel indices, 7/ and v are spin 
indices and is the 2x2 identity matrix. The trace in 
Eq. ([5a]) is taken over both sets of indices. We find that 
the average of the spin transmission vanishes in all cases, 



</d\ 



/RMT 



= 0. 



(9) 



For the (3 = 4 ensemble, this result was first obtained in 
Ref. Q • We further obtain 



varprj^^SRS 



= 4 



var[T) 



NiNjN^N^ 

N(4N 2 - 1) ' 
NiNj(N - 1) - NtN? 



(10a) 
(10b) 

(10c) 



N(2N- l)(27V-3) 

Eq. (fTUcj) first appeared in Ref. 0. We see that Eqs. Q, 
(|10b[) and (jlOcp all agree in the limit Ni j 3> 1, however, 
while the semiclassical expression Eq. (|7|) is valid only in 
that limit, Eqs. (|10[) are exact for any number of channels. 
Most interestingly, for a two-terminal setup with JVj = 
1, Eq. (Unci) gives var[T^0 =4 = and thus G fl = 0, 
in agreement with Ref. [16j . This restriction no longer 
applies once TRS is broken, as reflected in Eq. (|10bl) 
breaking TRS can turn spin currents on. 

Numerical simulations. We numerically confirm 
our findings using the spin kicked rotator model [13] • It 
is represented by a 2M x 2M Floquet matrix fl8j | 



2r w = (IK/XE/tllV, J, J' = 0,1, 



,,M 



M ' e cro, 



Xr, 



Swe 



-i(M/4w)V(2wl/M) 



1, (11a) 

(lib) 

(lie) 

(lid) 



The matrix II represents free ballistic motion, pe- 
riodically interrupted by spin-independent and spin- 
dependent kicks given by the matrix X, and correspond- 
ing to scattering at the boundaries of the quantum dot, 
as well as SOI. We choose 

V(p) = K cos(p + 9) cro + K so (a x sin 2p + a z sinp) . (12) 

The map is classically chaotic for kicking strength 
K > 7.5, and K so is related to the SO coupling 
time t so (in units of the stroboscopic period) through 
t so = 32tt 2 / 'K 2 M 2 From ([TT]). we construct the 

quasienergy-dependent scattering matrix as 

5(e) = P[e~ le - J(l - P T P)]~ 1 JP T , (13) 

with P a 2N x 2M projection matrix 

r 5 aP if/t' = /( fe ), 



P 







otherwise. 



(14) 



The l( k) (k = l,2,...,2iV, labels the modes) give the 
position in phase space of the attached leads. The mean 
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Figure 2: (Color online) Weak localization corrections to 
(top), and variance of (bottom) the charge (empty symbols) 
and spin (full symbols) conductance for the two-terminal 
kicked rotator of Eqs. (|11[) . Parameters are td = 10,20, 
K = 40, 60, 80, 90, K so = 120 K BOC and M = 128, 256, 512. 
The dashed lines indicate the RMT predicted crossover from 
fit — 4 to p = 2 [T3]. Our semiclassical prediction of Eq. (J7J) 
is illustrated by the straight black line in the bottom panel. 
For all data, TV > 10. 

dwell time td is given by td = M/N. The parame- 
ter K so breaks SRS over a scale K soc = 4iry/2/MT^/ 2 

and 9 breaks time-reversal symmetry over a scale 9 C = 
1/2 

4-7T / KMtjJ when Zo is finite. In our numerics we fix 
l = 0.14. When K > 1 and 9/9 c > 1, the charge 
conductance properties are those of the (5 = 2 ensem- 
ble, while for 9 = and K so /K soc 3> 1 they are those 
of the j3 = 4 ensemble [T3] - In °ur numerics, we fix 
K so /K soc = 120 and vary 9 to gradually break TRS, 
starting from 9 = 0. For simplicity, we specify to two- 
terminal setups and accordingly calculate the dimension- 
less spin conductance defined by Eq. ([3]) as = f° r 
(i = z. We checked, but do not show, that numerical re- 
sults remain the same if instead we consider [i = x,y. 

Fig. [5] first shows data for quantum corrections to the 
charge and spin conductance, as TRS is gradually broken. 
The top panel shows that weak localization corrections 
to the charge conductance are damped by a Lorentzian 
~ [1 + (919c) 2 }' 1 as predicted by RMT [| and semi- 
classics [lfj. There is no weak localization correction 
to the average spin conductance, both with and with- 
out TRS, in agreement with Ref. The bottom panel 
shows that charge conductance fluctuations are halved 
upon TRS, breaking and their behavior once again agrees 
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Figure 3: Spin conductance fluctuations for the kicked 
rotator with SOI defined in Eq. vs. the rescaled TRS 
breaking parameter 8/8 c for JVr = Al = 1. For 9 = 0, one is 
in the p = 4 ensemble and TRS forces the spin conductance 
to vanish [l(|. Breaking TRS results in a finite variance 
of the spin conductance. Dashed line: RMT prediction 
var[G M ] = 4/30 for Ar = Nl = 1 [see Eq. (flOE)) ]. Data 
correspond to K = 45, K so = 384(7r/Af) 2 with M = 128 
(red circles), 256 (blue triangles) and 512 (black diamonds). 
The curves do not lie on top of one another, because the 
rescaling of the horizontal axis with 8 C assumes Nb,,l 2> 1 QjJ- 



TRS 


SRS 


Charge transport 


Spin transport 


Yes 


Yes 


P = l 


P = 1; G M = 


Yes 


No 


= 4 


P = 4; Eqs. © and (fTOcl) 


No 


Yes 


P = 2 


G M =0 


No 


No 


p = 2 


Eqs. 1(9)1 and pObl 



Table I: Universality behavior of charge and spin transport 
properties in the four possible cases of broken or unbroken 
SRS and TRS. When both symmetries are broken, the spin 
transport properties correspond to those of the P = 4 Dyson 
ensemble in the limit Nn,Nh 3> 1. Deviations from P = 4 
are given in Eq. (|10p for the spin conductance variance. They 
are largest for small number of channels. 

is broken, see Eq. (jlObp . Spin and charge transport uni- 
versality classes are related to TRS and SRS in Table. Q] 
Examples of systems with broken SRS and TRS include 
spin-orbit coupled systems under external magnetic fields 
that we discussed, but also systems with spin textures 
and even spin valves with non-aligned magnetizations. 

We thank M. Biittiker for several interesting discus- 
sions at various stages of this project. 



well with theoretical predictions. The situation is entirely 
different, however, for the spin conductance fluctuations, 
which are essentially insensitive to the breaking of TRS. 
This is in agreement with our predictions, Eqs. ([7]) and 
(ITU|l for the large number of channels A^ > 10 considered 
in all data in Fig. [5] The new universal behavior corre- 
sponding to broken SRS and TRS emerges at larger 0, 
where the charge conductance corresponds to the /? = 2 
Dyson ensemble, while the spin conductance is very close 
to those of the j3 = 4 ensemble. 

Fig. [3] best illustrates the new universal behavior. 
When the exit lead carries a single transport channel, 
TRS requires that the spin conductance vanishes , re- 
gardless of the presence or absence of SRS. Fig. [3] shows 
that, when SRS is broken, breaking TRS turns spin cur- 
rents on, whose variance is given by Eq (|10b[) once TRS 
is totally broken. Note that the magnitude of the field 
necessary to break TRS for A/r^l = 1 becomes smaller 
and smaller in the semiclassical limit, M oo as the 
dwell time grows in that limit, td ~ M. 

Conclusions. By direct calculation we have shown 
that the spin conductance is an observable that is sensi- 
tive to the presence or absence of SRS even when TRS is 
broken. Breaking of SRS is necessary to magnetoelectri- 
cally generate a spin current, thus to acquire a finite spin 
conductance, but the latter is affected by TRS only when 
there are very few transport channels. Accordingly, we 
conclude that the j3 = 2 universality class splits into two 
different subsets for spin transport. In both cases, charge 
transport properties correspond to the /3 = 2 class, how- 
ever, the spin conductance vanishes identically when SRS 
is preserved, but exhibits a universal behavior when RS 
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